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Abstract 

We present a simple model for the doped compound Nd2-yCe y CuO±, in order 
to explain some recent experimental results on the latter. Within a Hartree-Fock 
context, we start from an impurity Anderson-like model and consider the mag- 
netic splitting of the Nd-Af ground state Kramers doublet due to exchange inter- 
actions with the ordered Cu moments. Our results are in very good agreement 
with the experimental data, yielding a Schottky anomaly peak for the specific 
heat that reduces its amplitude, broadens and shifts to lower temperatures, upon 
Ce doping. For overdoped compounds at low temperatures, the specific heat be- 
haves linearly and the magnetic susceptibility is constant. A smooth transition 
from this Fermi liquid like behavior ocurrs as temperature is increased and at 
high temperatures the susceptibility exhibits a Curie-like behavior. Finally, we 
discuss some improvements our model is amenable to incorporate. 
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We develop a model for the compound Ndi- y Ce y CuO^ following both the conjec- 
tures and the experimental results as given by There, a Schottky-like peak in the 
specific heat curve C(T) was observed. As Ce doping is increased, this peak reduces its 
amplitude, broadens and shifts to lower temperatures. In the overdoped case, at low 
temperatures, C(T) behaves linearly and the magnetic susceptibility remains constant. 
At high temperatures, the former behaves exponentially and the latter in a Curie-like 
fashion. A large effective mass was inferred and interpreted as a heavy fermion fea- 
ture. Czjzek el al. [ l] interpretation is as follows: the Schottky-like peaks result from 
the magnetic splitting of the Nd-4f ground state Kramers doublet due to exchange 
interactions with the ordered Cu moments. Taking at face value their interpretation 
and data, we develop a simple model for this compound's behavior. 

We start with an impurity Anderson-like model 0, within a Hartree-Fock {HF) 
context, and follow a somehow similar modeling as in ||. We denote by e k and n ku = 
c \ a c kcr the conduction band dispersion and occupation number operator, respectively. 
For the undoped case the Fermi level (taken as zero throughout this calculation) lies 
in the gap below the bottom of the conduction band (the undoped compound exhibits 
semiconducting behavior and undergoes a metal- insulator transition at y = 0.13 [|TJ). 
As the system is doped with electrons Q the band 'moves' towards the Fermi level (we 
assume the metal insulator transition of the band crossing type ||). In the present 
modeling, we only consider the doublet /-electron state structure, labeled by the spin 
variable a (orbital degeneracy and crystal field effects are neglected at the present time). 
Let rifi a = fj a fia be the localized /-electron level occupation number, at energy E Q . 
The / levels are assumed singly occupied (equivalent to consider the /-intra atomic 
Coulomb term Uff — > oo). Furthermore we denote by J the exchange interaction 
between the / and the conduction electrons. The Cu moment coupling via exchange 
to the /-electrons is modeled in a non-self consistent fashion, and the Cu moment mean 
value times the exchange constant is denoted by D. Let h be the external magnetic 
field, and the conduction (c) and /-electron magnetic moment mean values be denoted 
by Tn c j, respectively (absorbing all magnetic constants). Finally, the conduction-/ 
electron hybridization is modeled via the Anderson impurity model 0, so the HF 
Hamiltonian for this model is cast as 



except for a constant term and where V is the hybridization constant, N the total 
number of sites, and 
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Eka = £k — cr(Jm,f + h) E fa = E Q - a(D + Jm c + h) 



Notice the similarity of this Hamiltonian with the one presented by [[| in a related 
but somehow different context. Following the standard HF approach, as given by 0, 
we obtain the selfconsistent equation for the /-moment 



Imy-^ f du)f T (u)) 



7r ' J uo — Ef a + zA(u; 

with /t(w) the Fermi-Dirac distribution and the /-level broadening factor A is a func- 
tion of the energy u, the hybridization constant V and the conduction electron density 
of states as given by A(u) = 7rV 2 p c (u) [0. Notice the /-moment is normalized to 
the effective moment /Lt e // @- For a wide band metal, A is well approximated by its 
value at the Fermi level. Here, this is not the case since a metal insulator transition 
occurs as the compound is doped. We approximate A by its energy average (around 
the Fermi level), as A (a;) ~ A(y) = nV 2 p c (y). This 'average' is a monotonic func- 
tion of Ce doping. The conduction band moment can be likewise approximated as 
m c pa 2(Jrrif + h)p c (y). Once the / -moment is computed (selfconsistently), we are in 
a position to compute the / contribution to the energy as 



_Imy^ f duuf T (uj) 

tt u-E f(T (y)+iA(y) 1 } 

The specific heat and magnetic susceptibility are readily computed, at zero magnetic 
field, respectively as 



Our results are summarized in the figures. In Figure 1 and 2 we plot the spe- 
cific heat (in arbitrary units) and the specific heat over temperature, respectively; in 
Figure 3 the / magnetic susceptibility and in Figure 4 the / magnetic moment; all 
versus temperature in log scale. We attempted a quantitative fitting to the specific 
heat experimental curves ffl. In each Figure we plot four cases, corresponding to the 
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parameters' configurations (in Kelvin units) listed in the Table, with y the Ce doping 
relative concentration. The effective mass is estimated via the expression 



p f (e F = 0, T = 0) = ^1 j 

Ai(y) + (D(y) + J(y)m(y)) 2 



m* 2h 2 , 

— = -p f £F = 0,T = 

m mcr 

where a is the n.n. Nd lattice constant, and taken from [[I]]. In the last expression we 
fit the Lorentzian /-broadening with an effective tight binding Hamiltonian density of 
states. p[ . Furthermore, we consider E Q — e F ps (a non selfconsistent approximation, 
used in our model to approximately 'force' a singly occupied / -level). The parameters 
are computed within a virtual crystal approximation scheme f| |10| . This linear 
interpolation scheme mimics the effect Ce doping, thus A must increase with doping, 
as electron injection enhances the conduction density of states at the Fermi level. 
Correspondingly, electron doping must thin down the 'magnetic' parameters D and J, 
respectively associated to the incomplete filling of the Cu and Nd orbitals. 

From Figures 1 and 2 we notice the overall linear behavior of the specific heat at low 
temperatures, and an exponential behavior at high temperatures, in agreement with 
the experimental data M. Also, as doping (A) is increased, the Schottky-like peak 
moves to lower temperatures, reduces its amplitude and broadens, again in agreement 
with the experimental data. Notice that within our crude modeling for A, a non linear 
behavior at low temperatures is evident only for near null values of A, nevertheless 
our modeling stands as a good description for the metallic phase of this rare earth 
compound (this last remark is also valid for the computed susceptibility). A(u,y) is 
a function of both energy and doping, in fact proportional to the conduction band 
density of states p c (uj, y). For the undoped case (insulator) p F = p c {uj = e f , y) = 0. As 
doping is increased, p F becomes finite at the insulator-metal transition. Consider our 
expression for the specific heat (equations (|3p and (^)) with the energy independent 
approximation A(y). For the case T C A, the resulting Digamma functions can be 
expanded (Sommerfeld expansion) rendering the well known linear expression C = 7T. 
On the other hand, for the case A = and with e = — Ef^\, it can be readily 
shown that 
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T) (l + exp(^)) ; 
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the classical Schottky expression (see for example |TI]), with an exponential behavior 
at low temperatures and a sharp maximum at T « e. So as A varies away from zero 
we expect a transition from an exponential to a linear behavior. Thus, our single 
particle model accounts for the above mentioned transition as can be inferred from the 
experimental curves |I[ . The approximation A(u,y) ~ pf is a good approximation 
irrespective of the value of pp provided p c (uj,y) is a slowly varying function. For low 
doping this is not the case (an insulator-metal transition takes place), therefore if an 
energy independent approximation should be made for A, equating the latter to pp 
turns out to be a naive guess. Nevertheless qualitative behavior is obtained; namely 
the exponential to linear behavior in the low temperature region for the specific heat, 
and the maximum peak of the latter reduces its amplitude, broadens and shifts to lower 
temperatures as doping (A) is increased. We chose the effective (energy independent) 
A's to best fit the experimental data on the specific heat, specially with regard to the 
peaks' height and width. The low temperature regime for the underdoped cases exhibits 
a minute non linear behavior (see Figs.l and 2). Detailed analysis of our results shows 
a very weak exponential behavior 'creeping in' as doping is reduced. We believe this 
quantitative shortcoming at both low temperature and doping is unavoidable within 
the effective A approximation. This situation can be remedied by considering in earnest 
a 'rapidly' varying A(u),y). 

Figure 3 displays the moment suppression as doping and/or temperature are in- 
creased. Figure 4 displays the Pauli-like behavior of the system at low temperatures, 
and a Curie-like behavior at high temperatures. This feature agrees with the experi- 
mental data, for the metallic case y = 0.20 (m* ~ 1000m). 

From Figures 2 and 4, and only for case d we notice a smooth transition from a 
Fermi liquid behavior (with a large effective mass ~ 1000m) to a (two level) free-ion 
behavior at high temperatures. All other cases (m* ;$ 600m), and for which no data 
was presented in Q, exhibit a maximum at intermediate temperatures, suggesting an 
intermediate valence profile ||. 

The specific heat peak occurs at the same temperature where the moment is sup- 
pressed. Its amplitude is related in a monotonic fashion to the maximum moment 
value. This feature cannot be understood by just tuning the /-level splitting. Thus, 
the doping dependent hybridization mechanism seems to be an essential ingredient to 
explain the observed Schottky-like phenomena. 

Two distinct regimes seem to be present. The first one, hereafter denoted interme- 
diate valence regime (IVR, composed of the a, b and c cases, with m* ;$ 600), and the 
second, associated to case d (m* ~ 1000m) with no susceptibility peak, as seen from 
Figure 4. Following ]1] we identify it as belonging to a heavy fermion regime (HFR). 
Both regimes can be discriminated by their effective mass ratio and by their respec- 
tive susceptibility profile (Figure 4). There, the IVR shows a maximum in the region 
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separating the Pauli behavior (low T) from the Curie behavior (high T). In our model 
calculation such a maximum is smeared out in the HFR and the /-moment is further 
reduced as we cross over from the IVR to the HFR. (see Figure 3). Consistent with this 
picture, spin fluctuations are large at the Schottky peak temperature for the IVR (see 
Figure 4). In contrast, the HFR exhibits, as temperature is lowered , monotonically 
increasing spin fluctuations, needed for the /-moment reduction. 

We conclude that in our model the Cu moment exchange with the /-electrons is 
responsible for the appearance of the Schottky-like peaks. The peaks' structure are 
driven by doping via a simple Anderson mechanism. The latter also shapes the linear 
behavior of the specific heat and the Pauli-like susceptibility, at low temperatures. As 
doping is increased the system seems to cross from an intermediate valence to a heavy 
fermion regime. All these features are in good agreement with the experimental data 
presented in ref. 0. 

Our model stands as a simple and relevant analysis of this Schottky like effect, 
within a Hartree-Fock scheme. We believe that in order to explore further results of 
this model, such as the temperature scale separating the Fermi liquid-like behavior 
from the high temperature free-ion behavior, or to attempt further quantitative fit- 
ting with experimental results, we must first incorporate correlation effects in earnest, 
and lift some of the approximations used in our model. We briefly mention some im- 
provements our model is amenable to incorporate. We should go beyond HF, as in 
p|,P|,[P!0[],[P^|-P^] so the Uff —> oo limit is adequately treated, as for example with 
an effective medium approach ||,||-[|10|]. Within this scheme, both non-selfconsistent 
restrictions used in our model namely: E — Ep Pd and to consider the Cu moment 
as a 'classical' parameter (D), can be readily lifted and treated in earnest. Further- 
more, the approximation for A(u,y) as solely a doping dependent function can be 
relaxed. See for example the phenomenological scheme in []IS| and the renormalization 
group approach in |16| . These above mentioned improvements will be considered in a 
forthcoming paper. 
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Table 





V 


A [°K] 


D [°K] 


J [OR] 


m* x 1Q 3 

m 


a 


0.00 


0.90 


0.50 


3.50 


0.23 


b 


0.10 


1.20 


0.38 


2.45 


0.34 


c 


0.15 


1.35 


0.31 


1.93 


0.57 


d 


0.20 


1.50 


0.25 


1.40 


1.14 



Figure Captions 

Figure 1: Specific heat (arbitrary units) versus temperature. Four parameters config- 
urations in Kelvin units, labeled a,b,c and d, respectively (see Table 1). 

Figure 2: Specific heat (arbitrary units) over temperature, versus temperature. 
Labeling follows Fig. 1. 

Figure 3: /-magnetic moments versus temperature. Labeling follows Fig. 1. 

Figure 4: /-magnetic susceptibility versus temperature. Labeling follows Fig. 1 
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